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ABSTRACT 

From detailed spectral analysis of a large sample of low-redshift active galac- 
tic nuclei (AGNs) selected from the Sloan Digital Sky Survey, we demonstrate- 
statistically for the first time — that narrow optical Fe II emission lines, both per- 
mitted and forbidden, are prevalent in type 1 AGNs. Remarkably, these optical 
lines are completely absent in type 2 AGNs, across a wide luminosity range, from 
Seyfert 2 galaxies to type 2 quasars. We suggest that the narrow Fe II-emitting 
gas is confined to a disk-like geometry in the innermost regions of the narrow-line 
region on physical scales smaller than the obscuring torus. 
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1. Introduction 

Broad Fell multiplet emission are prominent features in the optical spectra of most 
type 1 active galactic nuclei (AGNs), but narrow optical Fell lines, either permitted or 
forbidden, are rarely seen. To our knowledge, narrow optical Fell emission has been reported 
in only three sources to date: I Zw 1 (Veron-Cetty et al. 2004), Mrk 110 (Veron-Cetty et 
al. 2007), and SDSS J1028+4500 (Wang et al. 2008). The apparent weakness of Fell 
emission in the narrow-line region (NLR) might be attributed to iron, a refractory element, 
condensing onto dust grains, which sublimate in the broad-line region (Laor & Draine 1993). 
This is especially true for low-ionization lines such as Fe II that arise from partially neutral 
portions of the NLR clouds (Ferguson et al. 1997). On the other hand, narrow forbidden 
Fell emission lines in the near-infrared (NIR) are observed commonly in both type 1 and 
type 2 AGNs (e.g., Simpson et al. 1996; Mouri et al. 2000; Rodriguez-Ardila et al. 2004, 
2005). According to these studies, NIR [Fell] emission most likely originates from AGN 
photoionization, although additional contributions from shock heating cannot be excluded. 
This raises the following questions: are narrow optical Fell emission lines really absent in 
AGNs? If so, why? 

It is possible that narrow Fe II lines are actually prevalent in the optical spectra of AGNs, 
yet have just been neglected or mistaken as a part of the co-existing broad-line emission or 
stellar features. Because of the complex atomic structure of the Fe + ion, there are so many 
transitions that individual Fell lines, even if narrow, are highly blended with each other, 
effectively mimicking much broader features. The situation is further exacerbated in type 1 
AGNs by the ubiquitous presence of broad Fe II multiplets as well as other emission lines. 
Without knowing a complete identification list of the optical transitions and their relative 
strengths, it is hard to discern them. Fortunately, Veron-Cetty et al. (2004) have recently 
identified and measured all the narrow Fell lines present in the optical spectrum of I Zw 1, 
a well-known narrow-line Seyfert 1 galaxy. Their analysis made it possible for us to explore 
systematically narrow optical Fe II emission lines in AGNs. 

In this Letter, we report the discovery of the prevalence of narrow Fell emission lines, 
both permitted and forbidden, in the optical spectra of type 1 AGNs and their nondetection 
in type 2 AGNs. We adopt a cosmology with H = 70 kms -1 Mpc" 1 , fi m = 0.3, and = 
0.7. 
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2. Sample and Data Analysis 

2.1. The Type 1 AGN Sample 

The type 1 AGN sample consists of the 4178 Seyfert 1 galaxies and quasars from Dong 
et al. (2010), selected from the spectral data set of the Sloan Digital Sky Survey Fourth 
Data Release (SDSS DR4; Adelman-McCarthy et al. 2006). Sample definition and data 
analysis methods are described in detail in that work. Here we provide a brief description 
and some additional treatments for verifying the existence of narrow-line Fell emission. 
Briefly, we select broad-line AGNs with prominent Fell emission in the SDSS bandpass and 
with continuum and emission lines suffering minimally from contamination by host galaxy 
starlight. The criteria are as follows: (a) z < 0.8, to ensure that Fe II emission in the restframe 
wavelength region 4434-4684 A (hereinafter FeIlA4570), as well as H/3 and [O III] AA4959, 
5007, are in the bandpass; (b) a median signal-to-noise ratio (S/N) > 10 per pixel in the 
optical Fell region (4400-5400 A); and (c) absence of detectable CaK (3934 A), CaH + He 
(3970 A), and H5 (4102 A) stellar absorption features at > 2a significance. 

At optical wavelengths numerous Fell multiplets and other broad emission lines are 
heavily blended to form a pseudocontinuum. Following Dong et al. (2008), we fit simultane- 
ously the AGN featureless continuum, represented by a power law, the Fe II multiplets, and 
other emission lines in the range of 4200-5600 A using a code based on the MPFIT package 
(Markwardt 2009). The Fell emission is modeled with two separate sets of templates in 
analytical formsp one for the broad-line system and the other for the narrow-line system, 
constructed from measurements of I Zw 1 by Veron-Cetty et al. (2004), as listed in their 
Tables Al and A2. Within each system, the respective set of Fell lines are assumed to have 
no relative velocity shifts and the same relative strengths as those in I Zw 1. In the optical 
spectrum of I Zw 1, the broad Fell lines are blueshifted by 150 km s _1 , whereas the narrow 
Fell lines are located at the systematic redshift of the host galaxy; this velocity separation 
made possible the identification of the narrow component from the broad one (P. Veron 
2009, private communication). The broad H/3 line is fitted with as many Gaussians as sta- 
tistically justified. All narrow emission lines are fitted with a single Gaussian, except for the 
[Olll] AA4959, 5007 doublet, each of which is modeled with two Gaussians, one accounting 
for the line core and the other for a possible blue wing, as seen in many objects. The redshift 
and width of narrow-line Fe II are set as free parameters. 

In order to verify that the narrow-line Fell are not residuals from poor broad Fell 



1 The implementation of the template functions in Interactive Data Language (IDL) is available at 
http://staff.ustc.edu.cn/~xbdong7Data_Release/FeII/Template/ . 
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subtraction due to mismatch of the broad-line Fell model, we use three schemes to model 
the profile of individual broad Fell lines: (A) the best-fit broad H/3 profile, (B) the best-fit 
total (narrow + broad) H/3 profile, and (C) a single Lorentzian with width set as a free 
parameter. Scheme A is a natural choice that has been used in our previous studies (e.g., 
Dong et al. 2008, 2010). Scheme B is rather extreme and formally unreasonable because 
H/3 definitely has a component from the canonical NLR; we use it as a stringent test of the 
reality of the narrow Fe II lines. Scheme C is inspired by the fact that the broad Fe II lines 
in I Zw 1 are best described by a Lorentzian profile (Veron-Cetty et al. 2004). The redshift 
of broad Fe II is set to be a free parameter. We note that the width of broad Fe II lines may 
be different from that of broad H/3 (see, e.g., Hu et al. 2008), yet we find that the flux of 
both broad and narrow Fell are quite insensitive to the exact profile assumed for the broad 
component (Dong et al. 2010). As noted by Vestergaard & Peterson (2005) and Landt et 
al. (2008), the broad Fell multiplets are so highly blended that the overall profile mainly 
depends on their relative strengths. As a strict demonstration that the detection of narrow 
Fell is robust, we show (§3.1) that we can recover these features in the residual spectrum 
even after excluding the narrow-line Fe II template from our model altogether. 



2.2. The Type 2 AGN Sample 

The type 2 AGN sample comprises the ~ 27, 000 Seyfert 2 galaxies in SDSS DR4 
selected according to the criteria of Kauffmann et al. (2003), having H/3, [Olll] A5007, Ha, 
and [N II] A6583 detected at > 5 a significance. The mean and standard deviation of their 
redshifts are 0.10 and 0.05, respectively. We subtracted the starlight continuum to obtain a 
clean emission-line spectrum using the method of Lu et al. (2006), using stellar templates 
broadened and shifted to match the stellar velocity dispersion of the galaxy. The stellar 
absorption lines must be subtracted well to ensure reliable measurement of weak emission 
lines. In the present study, we add the narrow Fe II template described above into the model, 
in order to detect any possible narrow Fell emission. Next, we fitted emission lines with 
Gaussians using the code described in detail in Dong et al. (2005); errors are given by 
MPFIT. There are 2671 objects in the galaxy catalog having a broad Ha component with 
S/N > 5; we excluded these from the type 2 AGN sample (cf. Zhang et al. 2008). 
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3. Results 

3.1. Narrow Fell Emission in Type 1 AGNs 

As a first step toward ascertaining whether narrow Fe II lines can be discerned without 
the aid of a narrow Fell template, we model the pseudo continuum with a power-law AGN 
continuum plus the broad Fell template only. For this purpose, we try all three schemes of 
modeling the profile of broad Fe II lines described in §2.1. The fitting results are very encour- 
aging: narrow Fe II lines are present in many residual spectra regardless of the broad Fe II 
model adopted. Figure 1 illustrates the three schemes applied to SDSS J113541. 19+002235.4. 
Plotted as vertical lines are the locations of the narrow Fell lines, both permitted (dotted) 
and forbidden (dashed), as well as some narrow transitions of CrII, Nill, and Till (cyan), 
that were identified by Veron-Cetty et al. (2004) and appear unambiguous in all the residual 
spectra. Those lines are also labeled in Figure 1. While there are small differences among 
the three fitting schemes, overall the pattern of narrow features appears robust. Moreover, 
the relative strengths of the narrow features seem, to first order, reproducible from object to 
object. This is shown in Figure 2 (top), where we plot the residual spectra for four sources 
fitted using Scheme A (broad Fell modeled using the broad H/3 profile). 

Next, we fit the entire type 1 sample with the narrow Fell template included into 
the pseudocontinuum model, as described in §2.1. The broad Fell line profile is again 
modeled with the above three schemes. The addition of the narrow-line component is justified 
statistically by the F— test p with over half of the spectra having a chance probability P nu \\ < 
0.01, regardless of the model adopted for the broad Fell profile (even the extreme case of 
Scheme B). Because the fluxes of both narrow and broad Fell are relatively insensitive to the 
choice of broad Fe II template, we adopt Scheme A (broad Fe II profile = broad H/3 profile) 
for the rest of the study. 

We detect narrow Fell emission at > 3 a significance in 2515 (60%) of the 4178 objects, 
and at > 5 a significance in 1872 (45%), using errors on the narrow-line Fell flux as given 
by MPFIT. This is consistent with the above F— test results. The restframe equivalent 
widths of the narrow Fe II features integrated over the region 4434-4684 A (Fe II N A4570) 
range from undetectable (^0.1 A) to 25 A. The intensity ratios of the narrow to broad 
FeIlA4570 feature, among the 2502 objects in which both components are detected at > 3 a 
significance, vary by 2 orders of magnitude, ranging from ~0.005 to 0.5, with a mean of 0.07 
and a standard deviation of 0.3 dex (computed in log-space). The equivalent widths of the 



2 We have found through experimentation that the F— test works well in such spectral fitting, although, 
theoretically, it holds only for linear models (see below and Dong et al. 2008). 
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two components do not correlate very tightly (Spearman correlation coefficient r s = 0.45; 
see Figure 3), suggesting that the narrow component is not an artifact of measurement 
uncertainty associated with deblending of the broad component. Remarkably, however, as 
described in the companion paper by Dong et al. (2010), the relative strengths of narrow Fe II 
with respect to the continuum and all other prominent emission lines in the near-ultraviolet 
and optical region (e.g., MgIlA2800, broad Fell, broad H/3, and [Olll] A5007) correlate most 
strongly with the Eddington ratic§| rather than with other physical parameters of the AGNs. 
A physical model unifying these and other correlations concerning Eddington ratio were 
proposed by Dong et al. (2009a, 2009b, 2010). 

It is difficult to obtain accurate measurements of the kinematics of the narrow Fell 
features. The lines are weak, most are severely blended, and both the spectral resolution 
and S/N of the SDSS data are limited. Notwithstanding these limitations, we make a 
preliminary statistical investigation for the 1032 objects with narrow Fell emission detected 
at > 10 er significance. The best-fitting narrow Fell model gives line widths (full width 
at half-maximum, corrected for instrumental broadening) mostly in the range of 200-800 
km s _1 , with a mean of 560 km s _1 and a standard deviation of 150 km s -1 ; these values are 
comparable to those of [OlIl]A5007 within ±50% (lcr). The velocity shifts of the narrow 
Fell system with respect to narrow H/3 range mostly from —200 km s _1 (blueshifted) to 
300 km s _1 (redshifted), with a mean of 50 km s _1 and a standard deviation of 100 km s -1 . 
With respect to the broad Fell system, the mean velocity shift of narrow Fell is 210 km s , 
with a standard deviation of 280 km s _1 . We cannot discern any obvious velocity differences 
between the permitted and forbidden narrow transitions. Without a robust estimation of 
the measurement uncertainties on these kinematic parameters at present, we refrain from a 
full discussion of the kinematics in this Letter. 



3.2. No Narrow Fell Emission in Type 2 AGNs 

In strong contrast to the case of type 1 AGNs, no Fe II emission is detected at > 3 a 
significance in any object in the type 2 sample! Considering that the detection of weak 
emission lines may be heavily influenced by uncertainties in starlight subtraction, we divide 
the type 2 sample into three subsamples according to [O III] A5007 luminosity, and build 
a composite spectrum for each. The composites are arithmetic mean spectra, constructed 



3 Eddington ratio, t = L/Z/Edd) is the ratio between the bolometric and Eddington luminosities. The 
Eddington luminosity (^Edd), by definition, is the luminosity at which the gravity of the central source 
acting on an electron-proton pair (i.e., fully ionized gas) is balanced by the radiation pressure due to 
electron Thomson scattering. 
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following Chen et al. (2010; see their §3.1). The luminosity ranges of the three subsamples 
are L[oni] < 10 4a5 , 10 40 ' 5 — 10 4L5 , and > 10 4L5 erg s _1 . To avoid possible confusion from 
the inclusion of star- forming galaxies, we further require the objects to be located above the 
maximum starburst line in the Baldwin-Phillips- Terlevich diagram (Baldwin et al. 1981) of 
Kewley et al. (2001), and to have [Olll] A5007/H/3 > 3 to avoid low-ionization sources. To 
obtain a comparison sample of type 2 sources with AGN luminosities more compatible with 
that of the type 1 sample (median L[oin] = 1-3 x 10 42 erg s _1 ), we also build a composite 
spectrum for the SDSS type 2 quasars presented in Reyes et al. (2008), which have 8 x 10 41 < 
L[om] < 4 x 10 43 erg s _1 . As shown in Figure 2 (bottom), it is clear that no narrow Fell 
features whatsoever are seen in any of the composite type 2 spectra. To quantify this striking 
result, we use the type 2 quasar composite, after starlight subtraction, to place upper limits 
on two relatively isolated features, FeIlA4925 (integrated in the vacuum wavelength range 
4918-4938 A, which is dominated by Fell 42 A4923 and FeIl]A4928) and Fell 49 A5234. 
Relative to [Olll] A5007, the limits are < 0.006 and < 0.005, respectively. By comparison, 
the mean values of these features in the type 1 sample are 0.05 and 0.02. 



4. Discussion and Summary 

The most surprising finding of this Letter is that narrow Fell emission is prevalent 
in type 1 AGNs, yet not present at all in type 2 AGNs. This is a statistically robust 
result, based on analysis of a large sample of individual spectra as well as very high-S/N 
stacked spectra. A possible explanation for this striking contrast, appealing to the canonical 
geometric unification scheme for type 1 and type 2 AGNs (e.g., Antonucci 1993), is that 
narrow Fe II emission arises from gas in the innermost regions of the NLR located interior 
to the obscuring torus. The inner edge of the torus has been estimated to be on scale of 
parsecs, roughly the dust sublimation radius (Suganuma et al. 2006); its radial extent is 
likely to be several tens of parsecs (e.g., Jaffe et al. 2004; see Granato & Danese 1994 for 
a model). In this picture the narrow Fe II-emitting region is visible along our line-of-sight 
in type 1 objects but obscured by the (extent of) the dusty torus in type 2 counterparts. 
As a low-ionization specie, Fe II may preferentially avoid the ionization cone and be largely 
confined to a disk-like geometry along the plane of the torus, as depicted in Figure 6 of 
Gaskell (2009). 

Yet, the complete absence of optical narrow Fell emission in type 2 sources appears 
perplexing in view of the fact that narrow [Fe II] emission in the NIR has been seen in AGNs 
of both types (see §1). Can line-of-sight dust obscuration be so effective in suppressing 
or hiding the optical lines, in all cases? Given the inherent patchiness of the interstellar 
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medium, this seems incredible. However, in the absence of a quantitative prediction for the 
intrinsic relative strengths of the optical and NIR Fell transitions, which, to our knowledge, 
has not been published, it is difficult to evaluate this hypothesis. Comparison between the 
NIR and optical lines may be further complicated by the fact that the transitions in these 
two wavelength ranges arise from different energy levels (< 1.5 eV for the NIR; ~ 3 eV for 
optical forbidden; ~ 5 eV for optical permitted), and hence likely trace physically different 
regions. The NIR lines probably tracer cooler, dustier material compared to the regions 
emitting the optical lines. 

There are several lines of fruitful work for the future. From an observational standpoint, 
it would be useful to obtain optical spectra of higher S/N and higher spectral resolution to 
better constrain the kinematics of the Fe II-emitting gas and its relation to regions emitting 
higher-ionization species. To better understand the puzzling mismatch between the optical 
and NIR lines, it would be instructive to obtain NIR spectra of the sources for which narrow 
Fell has been detected in the optical, and vice versa. And finally, theoretical calculations 
are needed to predict the intrinsic spectrum of narrow Fell across a wide wavelength range. 
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Fig. 1. — Fitting results of the optical spectrum of SDSS J113541. 19+002235.4, using three 
different models for the profile of broad Fell emission lines, without taking into account 
narrow Fell lines. We show the SDSS spectrum (black), the AGN power-law continuum 
(blue), the pseudocontinuum (AGN continuum plus broad-line Fell emission) with broad 
Fell lines modeled with the profile of broad H/3 (red), with the whole H(3 profile (purple), 
and with a single Lorentzian of variable width (green). In the bottom of every panel, the 
pseudocontinuum-subtracted residuals for the three respective models are displayed (with 
arbitrary offsets for clarify). Vertical lines and labels denote the permitted (dark gray, 
dotted) and forbidden (gray, dashed) narrow Fell lines, as well as some narrow CrII, Nill 
and Till lines (cyan), that appear unambiguous in the present spectrum. We also label 
several non-iron emission lines (black) and two broad Ti II lines (navy) . In the bottom panel 
we show the best-fit model for broad H/3 (gray, solid). 
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Fig. 2. — {Top) Residual spectra of four type 1 AGNs with strong narrow Fell lines, after 
the subtraction of the AGN continuum and the best-fit broad-line Fell emission modeled 
with the profile of broad H/3. The vertical lines are the same as in Figure 1. (Bottom) Four 
composite spectra of type 2 AGNs, ordered by [O III] luminosity. The narrow emission lines 
from non-iron elements are labeled according to the identifications by Vanden Berk et al. 
(2001) from their composite quasar spectrum. Strong emission lines are clipped for clarity. 
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Fig. 3. — Distribution of equivalent widths of narrow and broad Fell A4570 emission (black 
dots), along with 1 a error bars (gray), for the 2502 objects with both components detected 
at > 3 a significance. 



